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Abstract 

As three particles are advected by a turbulent flow, they separate from each other 
and develop non trivial geometries, which effectively reflect the structure of the tur- 
bulence. We investigate here the geometry, in a statistical sense, of three Lagrangian 
particles advected, in 2-dimensions, by Kinematic Simulation (KS). KS is a Lagrangian 
model of turbulent diffusion that makes no use of any delta correlation in time at any 
level. With this approach, situations with a very large range of inertial scales and 
varying persistence of spatial flow structure can be studied. 

We first show numerically that the model flow reproduces recent experimental results 
at low Reynolds numbers. The statistical properties of the shape distribution at much 
higher Reynolds number is then considered. Even at the highest available inertial range, 
of scale, corresponding to a ratio between large and small scales of L/rj « 17,000, we 
find that the radius of gyration of the three points does not precisely follow Richardson's 
prediction. The shapes of the triangles have a high probability to be elongated. The 
corresponding shape distribution is not found to be perfectly self similar, even for our 
highest ratio of inertial scales. We also discuss how the parameters of the synthetic 
flow, such as the exponent of the spectrum and the effect of the sweeping affect our 
results. Our results suggest that a non trivial distribution of shapes will be observed at 
high Reynolds numbers, although it may not be exactly self similar. Special attention 
is given to the effects of persistence of spatial flow structure. 

1 Introduction 

The transport of scalar fields by turbulent flows is an important process in many phys- 
ical situations ranging from the dynamics of the atmosphere and the ocean to chemical 
engineering. Specific examples of scalars are provided by pollutant density, temperature or 
humidity fields and the concentration of chemical and biological species pQ. 

Issues of transport and mixing in turbulence are directly related to the properties of fluid 
trajectories. The problem is thus often addressed using Lagrangian techniques [21 01 E] 
There is an established formal connection between the statistics of fluid particle motion 
and the concentration field of a diffusing scalar 0. Hence the important modeling issue of 
predicting passive scalar transport in turbulence can be addressed by following the evolution 
of Lagrangian particles [21 02 E] . 

The dispersion problem of one or two particles in the flow has been studied in great 
detail. In particular, the seminal work of Richardson |Hj leads to the prediction that the 
separation between two particles grows according to (R 2 ) oc et 3 , where e is the rate of 



energy dissipation in the now. Much less work has been devoted to the dispersion of 3 or 
more particles. The remarkable organization of the flow, which leads to the formation of 
very sharp fronts of scalar concentration |HJ EI] ; has a non trivial signature on the 
3-point correlation function of the flow [1111121 IT3*] . This, together with the well-established 
relation between the properties of the n-point correlation function and the properties of the 
evolution of n fluid particles advected by the flow, provides the motivation for studying the 
problem of dispersion of 3 particles or more. Despite recent progress both theoretically |14j 
and experimentally ^)JE3E1E]) little is known about the dispersion of 3 particles or 
more. 

An extra motivation to study dispersion of more than two particles comes from re- 
cent theoretical attempts to model turbulent velocity fluctuations in terms of small sets of 
Lagrangian particles ^S] 

The evolution of three particles configuration in turbulent flows has been considered 
numerically in direct numerical simulations (DNS) of 3-dimensional flows, at moderate 
Reynolds number [201 - Because of the limited range of inertial scales available in DNS, 
the numerical studies are unable to answer questions about shape statistics in the inertial 
range. A phenomenological model, introduced to describe the shape deformation in the 
inertial range, in high Reynolds number flows, lead to the prediction of a non trivial shape 
distribution [2U| . This model provided both the motivation and the theoretical background 
to analyze the experimental results of [2J- This experiment provided important insight on 
the statistics of deformation, although the experimental setup also suffered from the limited 
inertial range. 

In this paper, we consider the problem of Lagrangian dispersion of 3 particles with 
the help of the Kinematic Simulation (KS) method, introduced in [22| . KS provides a 
Lagrangian model of turbulent diffusion, based on a simplified incompressible velocity field, 
with a proper energy spectrum and no assumption of delta correlation in time made at any 
level. This model reproduces very well the Lagrangian properties observed in laboratory 
experiments [2H], as well as in DNS [21] ■ The computational simplicity of the KS allows us 
to consider very large inertial ranges : a ratio of scales of ~ 10 4 is easily accessible with 
moderate computer resources. KS thus turns out to be an ideal tool to study issues of 
dispersion in turbulent flows. 

This paper is therefore devoted to the study of Lagrangian dispersion of 3 points in a 
2 dimensional turbulent flow using the KS. We will demonstrate that the numerical results 
effectively reproduce the experimental results of |21j . and we will explore the large Reynolds 
number limit with the help of the KS. 

In section II, we briefly discuss the parameterization used to characterize the size and 
shape of the triangles, and we review the theoretical and experimental results on shape 
dynamics. Technical aspects of the simulation methods are described in section [21 The 
comparison between the experimental results of [21j and the KS simulations are presented 
in section 0] Section 03 contains our results concerning the large Reynolds number limit. In 
section we discuss the effects of persistence of flow structures on shape dynamics. Last, 
we present our concluding remarks in section VII. 

2 Previous results on shape dynamics 

In this section, we briefly review the previous work on shape dynamics |2()l 12 lj . 



2.1 Kinematics 

The evolution of a cluster of particles is described both by the overall scale, and by the 
shape of the object. In the case of a set of n = 3 particles, located at (i = 1, 2, 3), we 
define a set of reduced vectors involving relative separations only, defined by I12| : 
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The radius of gyration is defined as 
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where r?- = I Xj Xj I 9X6 the distances between the vertices of the triangle. The quantity 
R measures the spatial extent of the swarm of particles. In order to characterize the shape, 
we introduce a moment of inertia like tensor [2U] 

9 ab = jZptpl (4) 
i=l 

where p\ is the a-th spatial component of the vector p i . For a triangle in 2-dimensions, 
the tensor g has two eigenvalues, g\ > <?2 (note that 51 + 52 = R 2 )- These eigenvalues 
characterize the spatial extent of the swarm in the two principal directions. The ratio I2 
between the smallest eigenvalue, 52, and R 2 : 

* ^ § W 

provides us with a quantitative measure of the shape of the object. An equilateral triangle 
corresponds to I2 = 1/2. The smaller I2, the more elongated the triangle is. The moment 
of inertia tensor can be used both in 2 and 3 dimensions to characterize a set consisting of 
an arbitrary number of particles. 

In the case of a triangle, a full parametrization of the shape is provided by the quantities 
w and X) defined by : 
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w = 2 m_^ (6) 



By taking into account the symmetries of the triangle under any reparametrization of its 
vertices, the parameter space is restricted to < w < 1 and < x ^ 7r/6. The variables 
w and ^2 are related by the relation : I2 = (1/2)(1 — \f\ — w 2 ). A small value of w 
corresponds to a nearly collinear set of points. The quantity \ is small when the separation 
between two particles, say 1 and 2, is much smaller than their separation with the third 
one : r i2 << ri 3 ,r 2 3. 

We consider in this work the statistical properties of the shape distribution. To this 
end, we study the probability distribution functions of the various quantities R, I2, w and 
X characterizing the shape. The Gaussian distribution : Pc(pi,P2) = AAexp(— {p 2 + p 2 .)) 
provides an interesting distribution of reference. It can be shown (2^ that the distributions 
of x and w are uniform (in 2 dimensions) : Pg{x) = f and Pg{w) = 1. In particular, the 
corresponding mean values are : (x)g = J2, (f^G = \ and (12)0 = (1 — f )/2. 



2.2 Monte-Carlo model 



In order to study theoretically the distortion of sets of 3 or 4 particles by a turbulent flow 
in the inertial range of scales, a stochastic model based on phenomenological considerations 
was proposed in |2U| . At the heart of the model is a simplified scale decomposition of the 
full turbulent velocity field, on the scale of the global size of the triangle measured by the 
radius of gyration R [53]. Namely, the velocity field is written as: 

v = v<+v~+v>, (7) 

where v< is the contribution due to the small wavenumbers or large scales in the usual 
Fourier decomposition (\k\ < 1/2R), v> comes from the large wavenumbers (\k\ > 2/R) or 
small scales and v~ originates from the scales of the flow comparable to the global scale R 
(2/R < \k\ < 1/2R). The large scale contribution is uniform over the triangular configu- 
ration of particles, and is therefore assumed in [20] not to distort the set of particles. The 
v~ part of the velocity field acts coherently over the scale of the triangles with correlation 
time of the order of the characteristic time of turbulence at scale R, defined by : 

T (R) = i? 2 /3 e -l/3 ( 8 ) 

The small scale component v> is often assumed to be completely incoherent on the scale R 
of the three points and its correlation time is short compared to the characteristic time of 
turbulence at scale t(R). It is modeled in |20| by a white noise term. 

The action of v~ is approximated by a (coarse-grained) strain matrix, M ab = d a v b , 
acting on the vectors p^. The rapidly fluctuating, incoherent component v> is modeled by 
a Gaussian, white in time, random process. This leads to the following stochastic model 

M = p \M ab + ul (9) 
dM ab M ab 

~dT = ~^R) + ^ (10) 

where the indices i,j = 1,2 for 3 particles labels the relative vectors, see Eq. ©, and a, b 
labels the spatial components. The velocity fields u and the r] ab term are random Gaussian 
terms, delta-correlated in time with variances 

(Vab(t)Vcd(t')) = C 2 n 5(t - t') U ac 5 bd - ^5 ab sJ\ /t(R) (11) 

<tif(t)u}(0> = (y) 2 ^ " tr )^a b R 2 /r(R) (12) 

The stochastic model has been constructed in such a way that the matrix M is traceless 
(incompressibility) and correlated with a time scale t(R). Its amplitude is of the order of 
\M\ ~ ii -1 / 3 . The dimensionless parameter C v (respectively C v ) controls the importance 
of the incoherent jitter (respectively of the coherent term) in the model. 

Physically, the term p\M ab in equation (|10j) stretches and aligns the set of points. This 
distorting action is opposed by the action of the u term, which tends to make the shape 
distribution Gaussian. The shape distribution resulting from these two effects is non trivial, 
and depends continuously on the ratio C v /Cjj. A-priori, this number is of order 1. In the 
limit C v /C v — ► oo, the shape distribution becomes Gaussian. 



The model turns out to reproduce qualitatively several important aspects of the exper- 
imental results. A detailed analysis of the experimental data however pointed out several 
shortcomings of the stochastic approach |2J. In particular, it was found that whereas 
the stochastic model predicts a uniform distribution of the variable x, the experimental 
distribution of \ shows a peak near \ = 0- 



3 Kinematic simulations 

In contrast to the stochastic model described above, KS defines explicitly the velocity field 
which advects the particles. Following [2B], we define the turbulent velocity field v(x, t) by 
summing over a set of Fourier modes, k n : 

N k 

v(x, t) = ^ I A « A k " COs(k n • X + U n t) 

n=l 

+ B n A k n sin(k n • x + u> n t) (13) 

where is the number of modes in the simulations, k ra are the wave vectors, A n and B n 
are the amplitude vectors and uj n the frequency. The norms of the wavevector are chosen of 
the form k n = |k n | = kob n with a parameter b typically chosen to be b = (L/rf) 1 l <yNk ~ 1 \ The 
large (integral) scale, L, and the small (Kolmogorov) scale, r\ of the flow verify : L = l/k± 
and r\ = l/&/v fe {L/r] = b Nk ). The direction of the wave vector, k n = k n /|k ra | is uniformly 
distributed along the unit circle. Similarly, the directions of the vectors A n and B n are 
randomly distributed, and their amplitudes are chosen so that the energy spectrum is of 
the form E(k) ~ k~ p . The frequencies co n are taken to be uj n = \yJk'lE(k n ) : where A is a 
dimensionless parameter, a-priori of 0(1). The definition Eq. (|13|) makes the velocity field 
explicitly incompressible. Note that no delta correlation in time is used in KS at any level 
and that the parameter A controls the unsteadiness of the flow. 

To investigate the geometry of clusters of n = 3 Lagrangian particles, we simply advect 
numerically Lagrangian particles in the velocity field v(x, t) defined by Eq. ()13|) . This is 
done by solving a set of ordinary differential equations for the position vectors X(xo,i) : 

^X(xo,t)=v(x = X,t) (14) 
at 

with the initial condition X(xo,0) = xo. We start with an isotropic object, i.e. with 
an equilateral triangle, of a given size, and follow its evolution over time. The quantities 
characterizing the deformation of the object, such as I2, w, x are monitored as a function 
of time. We then perform ensemble averages over many triangles in different realizations of 
the velocity field to obtain the relevant particle statistics. 

We firstly validate the predictions from KS with the experimental results To this 
end, we choose the power of the spectrum p = 5/3, the ratio of inner to outer scales as 
suggested by the experiment, and we take A = 0(1). We then extrapolate our results to 
higher values of the ratio of inner to outer scales L/rj, to study dispersion in a high Reynolds 
number flow. 

We stress that KS is a model of the Eulerian velocity field, used to advect the particles. 
The KS velocity field, (|13|) has an interesting spatio-temporal flow structure, which varies 
with the parameters of the flow p and A [26,27]. Investigating systematically how the 
changes in the parameters p and A affects advection of particles is intrinsically interesting 
in the context of this study. 



4 Validation of KS 



Before making any predictions with KS regarding multi-particle statistics we firstly validate 
the model by comparing with the experimental results |21j . We are interested here in 
reproducing qualitatively the experimental results of |2*T] . This does not mean that KS is 
not able to reproduce quantitative predictions, at the cost of fitting the parameter A. We 
are merely interested in the trends of the shape evolution, that is, in the behaviors of the 
distributions of w, x etc as a function of time. 

An experimental investigation of the problem of dispersion of triangles by a turbulent 
flow was carried out in 2-dimensions, in the inverse cascade regime The flow was 

confined in a small container, 15 x 15cm 2 . Permanent magnets were placed under the 
bottom of the cell. The flow was stirred by running a current through a salted solution. 
The energy was injected at the scale Zj = 1.5cm. The velocity field was recorded by using 
standard Particle Image Velocimetry techniques, and was then stored on a 64 x 64 grid 
every 0.04s. The resulting spatial resolution was good enough to describe all the relevant 
scales of the flow. A Kolmogorov A;" 5 / 3 regime was observed, over the limited range of scales 
1.5cm < I < 5.5cm. The time resolution was also amply sufficient to follow numerically the 
evolution of particles. The evolution of a large number (~ 2 x 10 4 ) of triangles was then 
followed numerically. In this section we compare our results produced from KS with the 
experimental results of |21| and validate our model in the process. 

It was observed (see Fig. 3-5 of |21| ) that the typical size measured by the radius of 
gyration R of the triangles increased until it reached the largest scale 9cm of the ex- 
perimental setup where it started to fluctuate around this value. The evolution in time 
of the mean values of w and x ( see Figs. 4, 5 of (U) showed a rapid decrease of these 
parameters corresponding to strong shape distortions of the triangles. The smaller the 
initial separation ro the lower the minimal value of this parameter was observed. The 
shape distortion was maximum when R reached the lower value ^/ft? = 1.5cm of the 
inertial range. The mean values of these variables tend to an asymptotic value when 
R increases above the upper bound ^/~R? = 5.5cm of the inertial range. Specifically, 
it was found that (w) aS m = 0.5, {l2)asm ~ 0-11 and (x)asm ~ 0.26 These asymptotic 
values for w, I2 and x correspond to a Gaussian distribution of the p 1 and p 2 , which 
implies a uniform distribution for w, I2 and x with the Gaussian values {w)g<iu = 1/2, 
{h)Gau = (1 — 7r /4)/2 = 0.107 and {x)cau = vr/12 = 0.262 and a corresponding Gaussian 
distribution for R: P Gau (R) = (8R 3 /(R 2 ) 2 ) exp(-2R 2 / (R 2 )). 

The PDF of R and w, see Fig. 6 of 21 , can be well approximated by the Gaussian 
distribution for large values of time t = 80sec and t = lOOsec corresponding to values of 
the radius of gyration larger than the integral scale L. At later times, the finite size of the 
experimental system induces a saturation of the triangle size, so the tails of the distribution 
of R could no longer be correctly fitted by Pcau (R) EU ■ A very slow relaxation of the value 
of (x) towards its asymptotic, Gaussian value was observed. The fact that x = is more 
probable than x = 7r /6 implies that triangles with one edge much shorter than the two 
other ones has a large probability. This effect should ultimately disappear at later time, in 
the diffusive regime. 

The numerical experiment, consists of generating KS flows in two dimensions, with an 
energy spectrum E(k) ~ A; -5 / 3 , similar to the one observed experimentally, characterized by 
a ratio L/rj = 3.67, and with an unsteadiness factor A = 0.5. The smallest and the largest 



time scales of the flow are defined to be t v = 2tt/ Jk^E(k v ) and tg = L/u' respectively, 
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Figure 1: Time evolution of (R 2 ) 

where u' is the r.m.s velocity of the flow field. In this flow, we follow the evolution of 
three points, initialized as the vertices of an equilateral triangle of size r$. The results were 
averaged over ~ 10 4 configurations. 

The results obtained from KS show the same tendencies as the one observed in the 
experiment, as we now demonstrate. 

The behaviors as functions of time of the mean values of (R) (see Fig. (u>) (see 
Fig. |2I) and (x) (see Fig. EJ), as well as the evolution of the PDFs of R, w and % ( see Fig. 
4) computed from KS are very close to the ones observed experimentally (see Figs. 3-6 of 
|21j). Contrary to the laboratory experiment, where the spatial confinement of the setup 
induced a saturation of the growth of the radius of gyration, (R 2 ) grows like t at very long 
times. The distribution of sizes, R, as well as the distribution of w and x are Gaussian in 
this regime. The relaxation of the peak of probability for x ~ is faster in the KS than in 
the experiment. This small discrepancy is conceivably due to the large scale limitations of 
the flow. 

At intermediate time scales, when (i? 2 ) 1 / 2 is in the inertial range, the mean values and 
the distribution of w are very close to the ones observed in [2^ ■ 

These results demonstrate that KS reproduces very well the main properties of the evo- 
lution of three Lagrangian particles. A similar conclusion was reached in ,23 ; by comparing 
laboratory and KS results in two dimensions, and in [51] by comparing DNS and KS results 
in 3-dimensions. KS is thus a potentially very useful tool both for fundamental studies, and 
for dispersion studies in a more applied context. 

5 Predictions of KS in the large Reynolds number limit 

The results of the previous subsection demonstrate that the KS model reproduces quite 
satisfactorily the laboratory results concerning the evolution of three particles in a turbulent 
2-dimensional flow. We now investigate the large Reynolds number limit with the help of 
the KS. This is achieved by increasing the ratio of the largest to the smallest scale, L/rj 
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Figure 2: Time evolution of (w) 




Figure 3: Time evolution of (x) 



irrespective of initial conditions. 

A previous study has shown that the separation between two particles in the inertial 
regime shows sizable differences compared to the famous Richardson's scaling, according 
to which the separation grows according to (R 2 ) oc ei 3 irrespective of initial conditions. 
Specifically a strong dependence of two-particle dispersion statistics on the initial separation 
the particles was found [23] . I n the case of three particles studied here, it is of obvious 
interest to investigate the statistics of the radius of gyration of the set of particles. Also, 
the existence of a non trivial, time independent distribution of shapes, predicted by [2U] 
when the separation between particles is in the inertial range of scales, remains to be tested. 
We address these questions in turn, with the help of KS. 

We stress that there is no guarantee that the evolution of 3 particles is correctly predicted 
by the KS with p = 5/3 when the ratio L/rj becomes large. In spite of this uncertainty, the 
numerical results presented here are intrinsically interesting since the KS flow shares with 
real turbulence a number of important properties. 

The evolution of three particles has been followed numerically, using several KS flows, 
corresponding to different Reynolds numbers, or equivalently, to different values of L/rj 
(L/rj = 1691, 3381, 16909). Our runs, with the values of the parameters characterizing the 
flows, are listed in Table 1. 

The results of these simulations (Fig. 5, corresponding to L/rj = 1691 and Fig. 6, corre- 
sponding to L/rj = 16909) show some resemblance with the small L/rj case, see section 0J). 
As the Reynolds number is increased, the dependence of the variables describing geometry, 
(I2), (W), and (x) becomes weaker, both as a function of the initial size of the triangle and 
as a function of time. The values of these variables is always significantly different from the 
Gaussian values. 

The PDF of R, w and x & t the value L/rj = 1691 also show trends which are similar to 
the ones observed at smaller values of L/rj (Fig. 7). We have not followed particles long 
enough to see the Gaussian distribution of shapes at very large values of L/rj; the study has 
been restricted to the non trivial inertial range. In this range of scales, the variation of the 
PDF as a function of time is considerably weaker than observed in section and in this 
sense, the results suggest that one may be getting close to the self similar shape distribution 
predicted in 20 j. A visible deformation of the PDF of R can be observed as a function of 
time, reflecting the fact that at the later times, large excursions in the radius of gyration 
are getting close to the value of the largest available scale, L. 

Although still finite, the value of L/rj for run 13 is significantly larger than it is in any 
engineering, industrial or laboratory flow. Even so, our numerical results show persistent 
differences with the picture of a simple truly self similar regime. 

In all the reported DNS, in 2 [2Hj and in 3 dimensions |29j . a dependence of the behavior 
of (R(t)) on the initial separation tq has been reported. We observe a similar behavior in our 
KS calculation. Fig. 8 shows that the i 3 Richardson regime is never really reached. Instead, 
a continuous dependence of the variation of (R 2 )/t 3 on the initial separation ro is observed. 
This effect is clearly seen, even at our largest Reynolds number. This behavior suggests 
that the set of three particles always remembers its initial condition, which represents a 
departure with respect to the Richardson prediction. Also, at a fixed value of the initial 
size of the triangle, a non trivial power seems to emerge as the Reynolds number increases. 
A similar behavior was observed for the separation between two particles |23j . 

To investigate further this effect, we consider the PDFs of the radius of gyration, R. 
Fig. 9 shows a superposition of the PDFs of R/(R) for three different initial values of 
r o/v = 0-l> 1-0 5-0, and at different values of time. The distribution of the large values of 



Rj (R) are expected to be independent of time and Fig. 9 does not disprove this expectation. 
However, the distribution at small values of R/{R) seems to vary throughout the entire 
evolution. 

To compare different values of ro, we have plotted the PDFs of R/(R) corresponding to 
two different values of ro, but with similar values of {R) (Fig. 10, 11). For large enough 
times the PDF of R/{R) seem to collapse at large values of R/(R) within statistical errors. 
However, serious deviations are observed at small values of R/(R). In particular the peaks 
of the PDF at very small values of R/(R) are much sharper at smaller values than at larger 
values of ro. 

These results suggest that the distributions of R are influenced throughout the entire 
evolution by the initial value ro, insofar as a significant number of triplets do not really 
separate, and remain at a value R ~ ro. In the case of two particles this behavior has 
indeed been observed by Jullien et al. JH] m the laboratory for low Reynolds number flows, 
by Fung et al. j2H] and Nicolleau et al. [23] i n high Reynolds number KS simulations. In 
this way, the evolution depends in an essential way on the value of ro- 

The evolution of other geometric quantities, such as (w), (x) reflects to some extent the 
behavior described above (see Fig. 5, 6). Indeed, the lack of exact self similarity observed 
in the evolution of the radius of gyration, R, shows that the prediction of a truly time 
independent shape distribution is at best valid at Reynolds numbers impossible to attain. 
Although this prediction might constitute a good first order approximation, which becomes 
better as the Reynolds number increases, Fig. 6 shows that the mean values of w, I2 and \ do 
vary with time, even at L/77 = 16909. In addition, a systematic variation with ro is seen in 
Fig. 5, 6. The distributions are observed to remain non Gaussian as long as (R) remains in 
the inertial range, and do correspond to a higher probability of observing elongated objects, 
as anticipated in the stochastic model proposed in ^OHOl- In the light of the KS results, the 
stochastic model correctly predicts the main qualitative feature (the increased probability 
of elongated objects), but it doesn't incorporate the lack of self-similarity of multi-particle 
diffusion observed in KS. As we discuss in the next paragraph, memory effects relating to 
this lack of self-similarity are observed in laboratory experiments |1S| \21\. They are also 
observed in KS because, unlike stochastic models, KS incorporates the persistence of flow 
structures. Whether this lack of self-similarity persists in the laboratory and in nature at 
extremely high Reynolds numbers remains an open question. 

The results obtained so far suggest that the lack of self-similarity in the evolution of the 
radius of gyration, R, is due to the fact that particles stay together, at a distance of order 
ro with a high probability. Even if the radius of gyration grows to a value R ^> ro, it was 
observed in that two particles of the triplet can remain close to each other, with a large 
probability. We interpret the lack of a stationary distribution of shapes to reflect a similar 
cause : when one particle of the triplet separates from the two other ones, which remain at a 
mutual distance ~ ro, a very elongated shape is created. The relaxation of the distribution 
of shapes towards a stationary distribution will depend on how the two particles that are 
close together eventually separate. Our observations suggest that, the smaller ro, the longer 
particles will stay together, hence, the longer it will take the transient to relax. This lack 
of self-similarity is absent in stochastic models, such as the one proposed by ^F3\, and is 
consistent with the view that coherent streamline structures are persistent enough to cause 
a dependence of turbulent diffusion on initial conditions. 



6 Effects of persistence of the flow structure 



The KS model allows us to modify some of the characteristics of the advecting flow, both 
spatially and temporally. This is achieved by modifying the parameters A and p. The 
purpose of this section is to investigate the effect of changing the spatio-temporal structure 
of the flow, and in this way, to gain insight into the mechanisms involved in multi-particle 
dispersion. 

We firstly change the temporal structure of the flow by varying the persistence parameter 
A (see Fig. 12 ). This controls how fast the streamlines of the flow are jittered in comparison 
to the relevant eddy turn over time at the corresponding scale. This jittering makes the 
particles in the flow to be rapidly swept from one streamline to the other. Since in KS 
there is no interaction among modes of the velocity field, this jittering mimics the sweeping 
effects that are present in a real flow field. The minima of (hit)) and (w(t)) increase with 
increasing unsteadiness parameter A within the inertial range of time scales which means 
triangles are less elongated for larger values of A (Fig. 12). This happens because the paths 
of neighboring particles decorrelate faster for larger values of A and triangles quickly forget 
their memory of the initial state. Hence increased values of A should cause the triangle 
shape parameters to relax faster to their corresponding Gaussian values. 

Secondly we change the spatial structure of the flow field by changing the energy spec- 
trum i.e. changing the exponent p in E{k) ~ k~ p (see Fig. 13 ). This has the effect of 
changing the density of straining regions in the flow field |26| I27j thereby modifying the 
separation mechanism of particle pairs and clusters. The minimum of {^(t)) and (w(t)) 
decreases with increasing p (E(k) ~ k~ p ) within the inertial range of time and scales. This 
means that the clusters are more/less elongated during the inertial range of times when p is 
made larger /smaller. An explanation of this effect can be given in terms of randomness: as 
p increases there is less energy in the smaller scales of the turbulence which may mean less 
randomness leading to clusters remaining more elongated during the inertial range of times. 
In the context of KS as p increases there is indeed less energy in the small scales leading to 
smaller unsteadiness frequency uj n ~ kn p ^ 2 and therefore less randomness and more elon- 
gated clusters. However, a more searching explanation should invoke the spatio-temporal 
flow mechanism causing cluster elongation. One such mechanism already proposed in the 
literature |26l I22j is based on persistent effects of straining regions. The spatial density of 
straining regions decreases as p is made larger (2^1 E2j . 

However, the simple calculation, presented in the Appendix, indicates that the average 
straining rate per straining region increases faster than the number of such regions decreases 
when p is made larger. Provided that the effect of increased strain rate per strain region 
overwhelms that of the decreased number of such regions, then the same conclusion is 
reached: clusters should remain more/less elongated during the inertial range of times 
when p is made larger /smaller. 

7 Conclusions 

We have investigated the Lagrangian shape dynamics and the corresponding statistics of 
multiple particles namely of three particles advected by a two dimensional turbulent flow. 
We have used kinematic simulation (KS) to generate a turbulent velocity field and follow 
numerically sets of three particles in this flow field. The results of the simulation have 
been compared with the results of a two dimensional experiment We have identified 



a mechanism for the shape evolution of three particles depending on the underlying flow 
structure and the effect of persistence of these structures on the statistics of these shapes. 

Two regimes with well-characterized distributions have been identified in our simulation 
with KS. These regimes have also been observed in the two dimensional experiments of 
Castiglione et al. |21| . Two different regimes can be identified, according to the fact that 
the mean separation is large compared to the largest scale (diffusion regime) or in the 
inertial regime. The diffusion regime is characterized by a Gaussian shape distribution. 

In the inertial regime the scale rf -C {R 2 ) <C L 2 , we have observed what might be 
described as a trend, at best, towards a Richardson's law (R(t) 2 } oc t 3 . But the appearance 
of this regime was found to be dependent on the initial scale or size ro of the triangles. This 
has also been observed in the experiment |21j . 

The temporal evolution of (12(f)) , ( w (f) and (x(f) match with the experimental results 
of |2J. KS predicts the correct temporal evolution and the distribution functions of the 
above quantities. Monte Carlo simulation can do as well except for the PDF of x 

It is found that the clusters are more/less elongated during the inertial range of times 
when A is made smaller /larger (Fig. 12). The reason for this result must be that the 
persistence of the straining action of the flow is diminished when the flow is made more 
unsteady by increasing A. 

It is also found that the clusters are more/less elongated during the inertial range of 
times when p is made larger /smaller. An explanation of this effect can be given in terms of 
randomness: as p increases there is less energy in the smaller scales of the turbulence which 
may mean less randomness leading to clusters remaining more elongated during the inertial 
range of times. However, we also discuss a more searching explanation which invokes the 
straining mechanisms causing cluster elongation. 

The dependence of the shape of clusters on the initial separation between marked fluid 
elements is clearly demonstrated by the PDF of the radius of gyration of particle clusters 
with different initial separations not collapsing within the inertial range (see Figs. 10, 11). 
Our KS numerical experiments indicate that clusters tend to have memory of the initial 
state even when the turbulence has an extremely wide inertial range of more than four 
decades. 

This result is in agreement with the observed dependence |2.3j on the initial pair sepa- 
ration of the apparent power law governing the growth of inter-particle distances. If these 
effects are transient and due to a finite range of inertial scales, then our results indicate 
that they might only disappear at extremely high Reynolds numbers. However, the possi- 
bility should also be retained that these effects are not finite range transient effects, and 
are caused, instead, by the persistence of spatial flow structures at all scales, assuming this 
persistence remains at asymptotically high Reynolds numbers. 
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Appendix 



Assume we are given the spectrum E{k) = EoL(kL) p denned in the range 1/L < k < 1/rj 
of an isotropic turbulence. The mean square straining rate ((du/dx) 2 ) is proportional to 
Ii/l k 2 E(k) dk. Substituting the form of the spectrum and integrating we get for p < 3, 



In KS, the number density of straining stagnation points decreases with increasing p [2B] 
and [27j calculated the following scaling relation: 



where n s is the number density of straining stagnation points and D s is the fractal dimension 
of the spatial distribution of these points in the flow. In 2-D KS, D s = 3 — p |27| . Hence the 
number of straining stagnation points per unit area decreases with increasing p in our KS 
but the mean square strain rate per straining stagnation point, i.e. ((§f ) 2 )/ n si scales like 
£"0/(3 — p)L 2 . This implies that, although the number density decreases with increasing 
p, but the mean strain rate per straining stagnation point becomes stronger, which is the 
reason behind the decrease of the parameters (I2) and {W) with increasing p. 
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Figure 4: Time evolution of the PDF of R/a (a = \J (R 2 )), w and x( rac f° r r o/v = 1/6 in 2D 
KS, with L/r] = 3.67, /jb//;^ = 2.3 and A = 0.5. From top to bottom the figures are shown at times 
t = 2 x t, u t = 6 x 7j r) , t = 10 x t v and t = 14 x t n respectively. The light lines correspond to the 
Gaussian predictions P(i?) = 8(R/a) 3 exp[-2(i?/cr) 2 ], P(w) = 1 and P(x) = 6/tt [21]. 
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Figure 5: Time evolution of (R(t)), (I 2 {t)) 7 (w(t)) and (x{t)) produced by Kinematic simulation 
(KS) in 2D for L/rj = 1691 and A = 0.5 with t E /t n = 82.1. 
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Table 1: Different simulation parameters of our runs with kinematic simulation. 
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Figure 6: Time evolution of (R(t)), (hit)), (w(t)) and (x(t)) produced by Kinematic simulation 
(KS) in 2D for L/rj = 16909 and A = 0.5 with t E /t v = 538.1. 
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Figure 7: Time evolution of the PDF of R/a, w and x(rarf) for r /-q = 0.5 in 2D KS, simulation 
with L/rj = 1691 and A = 0.5 with tE/t v = 82.1. From top to bottom the figures are shown at times 
t = 20 x t v , t = 60 x t v , t = 100 x t v and t = 140 x t v respectively. 
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Figure 8: Time evolution of the (R(t) 2 ) / (t/t v ) 3 obtained by Kinematic Simulation of a triangular 
configuration of 3 particles in a two dimensional high Reynolds number (L/rj = 1691, 3381 and 16909) 
turbulent flow for different initial separation ro's Here the energy spectrum E(k) <~ fc~ 5 / 3 , A = 0.5, 
and the number of realizations is 5 x 10 3 . 




Figure 9: Time evolution of the PDF of R with different r 's obtained by KS in two dimensions 
with E(k) ~ fc^/ 3 , L/rj = 3381. From left to right r /?7 = 0.10, 1.0 and 5.0. 





Figure 10: PDF of radius of gyration R or global size with different values of ro- Here both ro and 
(R) have been normalized by rj. Two runs with different initial sizes are compared. The PDFs are 
shown at several times, for ro = 0.1 and r = 1.0. For comparison, the PDFs of R, corresponding 
to ro = 0.1 are shown at two times where the value of (R(t)} are the closest to the value of (R(t)) 
obtained with the larger value of r . The other parameters of the runs are L/rj — 3381, A = 0.5 and 
t E /t n = 130.2. 




Figure 11: PDF of radius of gyration R or global size with different values of ro- Here both tq and 
(R) have been normalized by rj. Two runs with different initial sizes are compared. The PDFs are 
shown at several times, for tq — 0.1 and tq — 1.0. For comparison, the PDFs of i?, corresponding 
to ro = 0.1 are shown at two times where the value of (R(t)) are the closest to the value of (R{t)) 
obtained with the larger value of ro. The other parameters of the runs are L/rj = 3381, A = 0.5 and 
ie/in = 130.2. The plots in this figure are the zoomed in version of figure ITU1 
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Figure 12: Time evolution of (hit)) and (w(t)) obtained by Kinematic Simulation of a triangular 
configuration of 3 particles in a two dimensional high Reynolds number (L/rj = 1691) turbulent flow 
for different A's. Here the energy spectrum E{k) ~ fc~ 5 / 3 , initial separation ro = 0.5 x rj and number 
of realizations is 5 x 10 3 . 
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Figure 13: Time evolution of (hit)} and (w(t)} obtained by Kinematic Simulation of a triangular 
configuration of 3 particles in a two dimensional high Reynolds number (L/rj — 1691) turbulent flow 
for different energy spectra E(k) <~ k~ p . Here A = 0.5, initial separation ro = 0.5 r\ and number of 
realizations is 5 x 10 3 . 



